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ABSTRACT 

We report on twenty-three clusters detected blindly as Sunyaev-Zel'dovich (SZ) decrements in a 
148 GHz, 455 square-degree map of the southern sky made with data from the Atacama Cosmology 
Telescope 2008 observing season. All SZ detections announced in this work have confirmed optical 
counterparts. Ten of the clusters are new discoveries. One newly discovered cluster, ACT-CL J0102- 
4915, with a redshift of 0.75 (photometric), has an SZ decrement comparable to the most massive 
systems at lower redshifts. Simulations of the cluster recovery method reproduce the sample purity 
measured by optical follow-up. In particular, for clusters detected with a signal-to- noise ratio greater 
than six, simulations are consistent with optical follow-up that demonstrated this subsample is 100% 
pure. The simulations further imply that the total sample is 80% complete for clusters with mass in 
excess of 6 x 10 14 solar masses referenced to the cluster volume characterized by five hundred times 
the critical density. The Compton y - X-ray luminosity mass comparison for the eleven best detected 
clusters visually agrees with both self-similar and non-adiabatic, simulation-derived scaling laws. 
Subject headings: Surveys - Cosmologyiobservations - Radio continuum: general - Galaxies: clusters: 
general - cosmic background radiation 
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1. INTRODUCTION 

Measuring the redshift evolution of the mass function 
of galaxy clusters provides powerful constraints on cos- 
mological parameters, complementary to those obtained 
from different measurements, such as the angular power 
spectrum of the co smic microwave background (e.g., 
iKomatsu et al.ll2010ft or luminosity d istances of Type la 
supernovae (e.g.. iHicken et al1l2QQ9[ ). The potential of 
cluster surveys as a cosmological probe has recently been 
demonstrated through the analysis of large X-ray and 
optically selecte d cluster samples spanning a wide range 
in redshift (e.g.JVikhlinin et al.ll2Q0llMantz etld1l2010l : 



IRozo et alJEjOlffl T 

A powerful method for cluster detection and charac- 
terization is provided by the thermal Sunyaev-Zel'dovich 
(SZ) effect. The SZ effect is the process by which en- 
ergy is transferred from the hot intra-cluster medium 
(ICM) of galaxy clusters to the cosmic microwave 
back ground (CMB) throug h inve rse Compton scatter- 
in gjZerdovich Sunvaevl 119691 : iSunvaev ZePdovichl 
1 19701 ). The overall effect is to distort the thermal 
spectrum of the CMB. In the case of the present 
study at 148 GHz, the SZ distortion produces a roughly 
arcminute-scale temperature decrement along the line 
of sight to a galaxy cluster. The SZ distortion is pro- 
portional to the pressure of the ICM integrated along 
the line of sight, providing a unique probe of clus- 
ter physics. The integrated SZ flux therefore scales as 
the thermal energy content of the cluster, which is ex- 
pected to scale with cluster mass for cluster gas hydro- 
statically supported by thermal pressure. The ampli- 
tude of the SZ distortion is approximately independent 
of the cluster redshift, in marked contrast to the ob- 
served cluster brightness in other wave bands. There- 
fore, the SZ effect is particularly suited for inventory- 
ing the distant universe. Galaxy cluster surveys are 
one of the few methods that probe the growth of struc- 
ture and interpretation of survey yields has the po- 
tential to precisely constrain the equation of state of 
dark energy (e.g.. iBartlett fc Silkl 119941: [Holder et al.l 
120001 : iHaiman et al.ll200lHMajumdar Mohrll2004fl . For 
a full description of the SZ effect, historical perspec- 
tives, and applications to cosmology , see the review art i- 
cles by Bernstein & Dodelson (1990), Birkinshaw (1991), 
iReph aeli (1995), Birkinshaw (1999), and lCarlstrom et al.l 
(|2002l ). 

The history of successful targeted observations of 
the SZ effect in previously known galaxy clusters goes 
back two decades. Pioneering measurements were mad e 
with the OVRO 40m telescope (Birkin shaw et al.lfl99lD . 
the Ryle Telescope pones et a J [l993|), the OVR O 5m 
telescope (jHerbig et al.l 119951 : iMvers et al.l 119971). th e 
OVRO/BIMA interferometers (jCarlstrom et al.l 119961 ). 
the bolo metric SUZIE re ceiver on the Caltec h Sub- 
millimeter Observatory (Holzap fel et al.l Il997| ). and 
the Nobeyama Rad io Observatory 45 m telescope 
(jKomatsu et al.| [l999). Scientifically these original stud- 
ies put special emphasis on deriving the Hubble param- 
eter from angular diameter distances to clusters mea- 
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sured through a combination of SZ and X-ray data. 
In the last decade, these geometric measurements have 
been expanded to statistically significant samples of 
clust ers (e .g., the OVRO/BIMA s tudies bv iReese et al.l 
120021 and iBonamente et al.l [2006) and serve as an in- 
dependent check on estimates that rely on the conven- 
tiona l extragalactic dista nce ladder (e.g. JFreedman et al.l 
120011 : iRiess et~a"T1 l2QQ9h . Beyond distance scale stud- 
ies, targeted SZ work has focused on understanding 
the SZ correlation to other cluster ob servables, es- 
pecia ll y those which tr a ce m as s fe.g.. iBensqn et aLl 



200 41: IBonamente et al.l 120081: iMcIrmes et all 200 



Marrone et all 120091: iMe lin et al .] 120101: [Huang et all 



2M0I : IPlagge et al J 12015 iCulverhouse et al.ll201of j. The 

SZ/X-ray and SZ/optical lensing correlations provide the 
mass calibration necessary for understanding the clus- 
ter mass function through cosmic time. By measur- 
ing the redshift evolution of the cluster mass function, 
constraints can be put on the matter density field and 
dark energy. Thus proper mass calibrati on is essential 
for st udying cosmology with clusters (e.g., Fr ancis et al.l 
|2005[ ). Finally, targeted SZ observations are expanding 
our knowledge of the physics of the I CM by probing the 
clust e r gas pressure at large radius (iMroczkowski et al.l 
120091 : iPlagge et all 120101: Komats u et al.l 120101) as well 



as at high resolution (Komatsu et al. 2001; Ma son et al.l 

HoTqD. 

The first blind detections of galaxy clusters through 
their Sunyaev-Ze r dovich effect were reported by 
iStaniszewski et al.l (12009ft fo r the South Pole Tele- 
scope (SPT; Carlstrom et al. 2009). Three of the clus- 
ters from Staniszewski et al. (2009) were previously un- 



known. iHincks et al.l (|2009f ) reported on the cluster de- 
tections for the A t acama Cosmology Telescope (ACT; 
iFowler et all 120071: l^tzejal.l I201O ). While the clus- 
ters from IHincks et~aTl ([2009) were previously known, 
they were blindly detected in the survey data. SPT 
and ACT have been able to realize blind detections from 
cosmologically significant surveys because they com- 
bine three essential design features: arcminute resolu- 
tion matched to the size of clusters, degree-scale fields 
of view for efficient surveying, and the unprecedented 
se nsitivity of 1000-el e ment bolometric detector arrays. 
In IVanderlinde et al.l (|2010l ). the SPT collaboration re- 
ported on twenty-one clus ter candidates, tw enty of which 
were optic ally confirmed (iHigh et al.l l2010"). In a follow- 
up study, lAndersson et al.l (j2010h presented the X-ray 
properties of the SPT sa mple and derived SZ/X- Ray 
correlations. Furthermore, IVanderlinde et all (|2010f ) ob- 
tained cosmological constraints from an SZ-selected sam- 
ple of clusters. 

In this work, the ACT collaboration reports on 
twenty-three clusters detected in a subset of the data 
taken during the 2008 observing season. This pa- 
per describes the detection and millimeter properties 
of the clusters. Companion papers describe the physi- 
cal properties and obser ved purity of the cluster sample 
(jMenante au et al. 2 O il]) as well as cosmological impli- 
cations (jSehgal et al.ll2010aT ). These papers and others 
cited above are part of a larger initial science release from 
ACT which also includes results on the CM B power spec- 
trum and associated parameter constraints ( Fowler et al 



2010; iDas et al.l 1201(1 iDunklev et aT1 l2010: Haiian et al 



2010[ ) and on compact millimeter sources (jMarriage et al 
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Fig. 1. — Sensitivity map with detections. The map shows the sensitivity over the subset of ACT 2008 148 GHz dataset considered for 
this study. The gray-scale encodes the noise rms in |oK of a map match- filtered for a /3-model (/5 = 0.86, C = 1.00'). The median noise in 
the map is 36 \xK. Black boxes mark the location of the twenty-three optically confirmed clusters. The size of each box is proportional to 
the corresponding cluster decrement. 



[20l0h . 

This paper is organized as follows. The data reduction 
and catalog construction are described in Section [2j The 
detections and their properties are discussed in Section [3l 
An estimate of purity and completeness of the sample is 
derived from simulations in Section HI Finally, in Section 
[5j the SZ signal is compared to X-ray-derived mass for a 
high significance subset of the cluster sample. 

2. FROM OBSERVATIONS TO CLUSTER DETECTION 

The observations and methods to reduce the raw data 
to maps for this work ar e the same as used in the 2008 
power spectrum study dFowler et all 201 Of) and extra- 
galactic source study ([Marriage et al.T2010[ ). The pri- 
mary difference between the source and cluster study 
lies in the compact signal profile used in the implemen- 
tation of the matched filter-based detection. Therefore, 
we restrict our discussion here to a short summary of the 
observations and d ata r eduction, and refer th e reader to 
iFowler et all (j2010l ) and lMarriage et al.l (|2010D for a more 
complete description. We describe in detail the particu- 
lar implementation of the matched filter for clusters. For 
a d escription of the A CT receiver and instrumentation 
see iSwetz et all (|2010f ). 

2.1. Observations and Data Reduction 

ACT is a six- meter telescope operating at 5200 m in 
the Atacama Desert of northern Chile. The site was cho- 
sen for its excellent atmospheric transparency and access 
to both southern and northern skies. The telescope has 
three 1024-element arrays of transition edge sensors op- 
erating at 148 GHz, 218 GHz, and 277 GHz. This study 
uses 148 GHz data from a 455 square-degree subregion of 
the 2008 southern survey. The subregion lies between 
right ascensions 00 h 12 m and 07 h 08 m and declinations 
— 56 !!' and — 49°00 / . Figure [T] is a map of the sensi- 
tivity across the subregion along with the locations of 
the clusters reported in this study. The median rms of 
the map optimally filtered for detecting a /3-model profile 
(/3 = 0.86, 6> c = 1.000 is 36^K. 

The data for this study were calibrated to the tem- 
perature of Uranus with a precision of 6%. The abso- 
lute positional uncertainty in the maps is established 



at 3.5" by compariso n of ACT radio source locations 
(M arriage "eTaj1l20T0l ) to cross-identified sources in the 
Australia Telescope 20 GHz Survey (jMurphv et al.ll201Q[ ). 
Note, however, that reported cluster locations have an 
additional uncertainty due to the effect of noise on the 
estimate of the cluster center. Considering the posi- 
tional uncertainty for sources detected be tween signal-to- 
noise ratio (SNR) 5.5 and 10 reported in iMarriage et all 
(|2010f ). we expect estimates of the ACT cluster center 
positions to scatter with an rms of roughly 10". Ad- 
ditional ambiguity in the cluster position arises in the 
case of an extended, non-spherical (e.g., disturbed) sys- 
tem. In such cases, the positional uncertainty can rise to 
arcminute scales. 

With calibration and astrometry solved, the final step 
in the data reduction is map making. We iteratively 
solve for the maximum likelihood (ML) estimate of the 
map using a custom precon ditioned conjugate gradient 
code (e.g. JPress et al.lll992f ). Because ACT samples the 
sky along multiple scan directions (i.e., the data is cross- 
linked), we are able to produce an unbiased ML map 
of the microwave sky with a faithful representation of 
structure from degree to arcminute scales. 

2.2. Filtering and Cluster Extraction 

For cluster detec tion we use matched filters 

(lHaehnelt fc Tegmarkl 119961 : iHerranz et all 120021 : 
iMelin et"aTI l2006f T We model the sky temperature 
fluctuation at a point x as 

ST{x) = ^ ST 0ii Be Cii (x - au) + ST other (x) (1) 

i 

where STo^ and Bq c ^ are the amplitude and unit- 
normalized profile of the i th cluster. In what follows we 
choose Bq c to be the isothermal /3-model ( j3 = 0.86 ) with 
core radius C ranging from 0.25' to 4.0' a nd convolved 
with a n isotropic ACT 148 GHz beam from lHincks et all 
(|2QQ9h . The profile is tapered to zero by multiplication 
with a cosine in the range 0.5-5.5# c . The choice of j3 
is motivated by the best fit to an average SZ profile in 
iPlagge et al.l (|2010f ) who found consistency between the 
/3-model fit and the Generalized NFW profile fit to Yx in 
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Fig. 2. — 148 GHz sub-map. The data have been weighted by a smooth function yj N t, s /N t, Sj max, where N t, s is the number of data 
per pixel. This weighting levels the amplitude of white noise across the map to that corresponding to the deepest data. The data is then 
match-filtered with a /3-profile (0 C = 1.0'). The coverage becomes shallower toward the bottom-left of the map causing a visible increase in 
the rms. The inset shows the flux density distribution across the data weighted by the square-root of the number of data per pixel. The 
data distribution is shown as a gray histogram on which is plotted a dashed Gaussian distribution with standard deviation 24|iK. Sources 
were removed prior to filtering. Three clusters are recognizable from left to right: the Bullet Cluster (bottom-left), A3404 (middle), and 
AS0592 (middle-right). 



lArnaud et all (|2009f ). The temperature field ST other (x) 
consists of noise modeled from difference maps, pri- 
mordial CMB flu ctuations with power spectrum from 
iNolta et"aLl ([2009f h as well as sources and undetectable 
(i.e., low mass) clusters. The spectral signature of the 
source and cluster contributi on to 5T n the,r ( x) is mod- 
eled from fits t o ACT data in iFowler et all (j2Q10l ). See 
iMarriage et al.l (|2010f ) for a full discussion of ST ot \ ier (x). 

Before filtering, bright (SNR > 5) sources are in- 
painted with sky temperature in the neighborhood of the 
source. Furthermore, the map is weighted by the inverse 
square root of the number of observations. This has the 
effect of flattening the white noise across the map. The 
map is filtered in the Fourier domain using a matched 
filter: 

S Bl{k')\ST M „{k')\-' B,,(k')Ak' 
where Bo c (k) and ST ot her(k) are the Fourier transforms 



of Bq c and ^T Q ther, respectively. The map is filtered us- 
ing /3-models with core radii from 0.25' to 4.0' in 0.25' 
steps. These core radii were chosen to span the range of 
angular scales characterizing massive clusters from low 
to high redshift. In a given map, the SNR of a detection 
is defined as the ratio of the extremum of the cluster 
decrement to the rms of the filtered map. The reported 
SNR for a given detection is the maximum SNR from the 
set of filtered maps. Figure [2] shows a subsection of the 
source-masked and filtered map containing three known 
clusters. Shown in an inset of this figure is the pixel 
flux distribution of the filtered map. The flux distribu- 
tion suggests the filtered map is dominated by Gaussian 
noise. 

The filter in Eq. [2] is normalized such that the central 
decrement in the filtered map is an estimate of the cen- 
tral decrement value in \±K. It is this value, recovered 
from the filtered map, w hich we use as our d ecrement 
estimate. As discussed in IMarriage et al.l (|2010f ). this es- 
timate suffers from increased error and bias due to the 
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Fig. 3. — Temperature decrements at 148 GHz of twenty-three clusters. Each thumbnail is approximately 8 ; on a side and is excised 
from the source-removed, filtered map. The grayscale of each thumbnail is symmetric about zero with a range set to the depth of the 
decrement. The ordering of the clusters is according to decreasing SNR from left to right and top to bottom. The structure northeast 
of ACT-CL J0641-4949 may be part of the cluster or may be contamination: the decrement of this nearby cluster may be noise boosted. 
ACT designations are at the top of each thumbnail and alternate names at the bottom, denoted "New" if newly discovered by ACT. 



0.5' pixel size of the ACT maps. In order to mitigate 
these errors, the pixel size of sub-maps around each de- 
tection is decreased by Fourier transforming the map and 
adding extra zeros to the boundary of the Fourier trans- 
form. The inverse of this zero-padded transform is then 
the original submap with smaller pixels. For this study, 
the pixel size is changed to 0.0625' which reduces the er- 
rors associated with the discrete approximation of Eq. [2] 
to less than 1%. The errors we quote for the map-derived 
decrements are simply the decrement amplitude divided 
by the SNR. 

3. CLUSTER DETECTIONS AND PROPERTIES 

Table [T] lists clusters detected by ACT. All detections 
in the table have been confirmed through optical follow- 
up imaging of a larger catalog of all SZ candidates from 
the ACT data (jMenanteau et al.ll2010h . The table is di- 
vided into two groups of candidates: candidates detected 
at SNR > 6 and candidates detected at lesser signifi- 
cance. The optical follow-up showed that for candidates 
detected at SNR > 6, the larger catalog of all SZ can- 
didates is 100% pure (See also Section HJ). Of the nine 
detections at SNR > 6, three are newly discovered sys- 
tems. Among the fourteen candidates confirmed at lower 
SNR, seven are new detections. Thumbnails of the clus- 
ter decrements are shown in Figure 03 



The multifrequency properties of the clu sters covered 
in thi s wo rk are discussed at len gth in iHincks et al.l 
(|2QQ9h and iMenanteau et all (|2010l ) . Paraphrasing these 
studies, X-ray luminosities were computed from Chan- 
dra and XMM observations where available and supple- 
mented by ROSAT All Sky Survey for the remaining 
clusters. These luminosities suggest that the full sample 
is composed of high mass clusters with a characteristic 
mass of M500C ~6 x 1O 14 M , where M500C is defined as 
the mass of the cluster at five hundred times the criti- 
cal density of the universe at the cluster redshift. The 
sample has a median redshift of 0.44. All but one of the 
clusters at z > 0.35 in the sample were discovered via 
the SZ. 

Of the new ACT clusters, the most notable is ACT- 
CL J0102-4915. This cluster has a decrement amplitude 
(ST = -1046±116/iiT) comparable to that of the Bullet 
Cluster (ST f = — 734±64/iif ) as well as a smaller angular 
core radius, consistent with its higher redshift (%>hoto = 
0.75). The X-ray luminosity (0.1-2.4 keV) derived from 
ROSAT observations is 18.39 ± 6.49 x 10 44 erg s _1 which 
is also comparable to that of the Bul let as well as the 
other massive clusters in our sample ( Menanteau et all 
l2010h . The cluster richness is also consistent with that 
of a massive cluster (N ga i w 50). The strength of this 
decrement may result from merger activity, elongation 
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along the line of sight, or other decrement-enhancing ef- 
fects. Planned spectroscopic and IR follow-up observa- 
tions will determine the true nature of the decrement. 

Table Q] contains estimates of the cluster decrement am- 
plitudes STq based upon the matched filter described in 
Section E2J Also listed is the core radius C correspond- 
ing to the filter that exhibits the largest SNR for each 
cluster. Figure 2] shows the temperature decrement and 
SNR of two newly discovered clusters plotted as a func- 
tion of filter scale (i.e., the core radius). For some of the 
clusters, such as ACT-CL J0346-5438 shown in the fig- 
ure, the SNR is a relatively broad function of core radius. 
The corresponding uncertainty in the size of the cluster 
results from the blending of the cluster decrement with 
the background (e.g., CMB). Th is confusion f or sin gle- 
frequen cy data was obser v ed by iMelin et al.l (j2006[ ) as 
well as IVanderlinde et all (|2010f ). For clusters with a 
more compact and pronounced core, like ACT-CL J0102- 
4915, the cluster-background confusion is smaller and so 
the optimal filter scale is easier to determine from the 
data. Confusion with the background inflates C while 
contamination by detector white noise causes a cluster 
to appear more compact than it actually is. The latter 
effect can be explained by the fact that most of the power 
in the white noise field is at the smallest scales. Simu- 
lations (See Section 2J) show contamination for clus- 
ters detected outside the range 6 C = 0.75' — 1.75'. For 
this reason, C estimates reported below 0.75' or above 
1.75' are suspect. Incorporation of ACT's multiple fre- 
quencies will help distinguish clusters from background. 
Additional years of data will help reduce the fraction of 
white-noise-boosted detections. 

The STo recovered from the filter is the amplitude of 
the cluster decrement in the ACT map. As such, it is the 
amplitude of the intrinsic cluster profile convolved with 
the ACT beam. In Table [T] we also include an estimate 
of the decrement 5Tq = 5TqC~ 1 (0 c ) where we have cor- 
rected the raw STo by the smearing effect of the beam on 
the height of the corresponding /3-profile. This correc- 
tion for a /3-profile with C = V is 1.67. The correction 
is larger for smaller core radii. To estimate the uncer- 



tainty in the recovered 5Tq, the cluster extraction algo- 
rithm was run on a simulation with realistic ACT noise, 
CMB and clusters with sym metric /3-mode l profiles . The 
more realistic simulations of iSehgal et al.l (|2010bl ) were 
not used as they were not provided an estimate of the 
decrement amplitude. As the /3-model profiles in this 
simulation are matched to our filters, the error estimate 
from these simulations should be taken as a lower limit. 
With those caveats, the mean error on STq from the sim- 
ulations is 20% greater than the statistical error inferred 
from the SNR (i.e., systematic errors are 40% of the total 
error budget). This modest increase in the error budget 
was accompanied by a large scatter in recovered scale 
(^30% with respect to the input scale value) with a non- 
Gaussian tail towards larger recovered radii. This implies 
that while the physical scale of the cluster may be con- 
fused with the background the amplitude of the peak 
is relatively immune to this confusion. When the cluster 
extraction algorithm was run on the /3-model simulations 
with prior knowledge of the scale, the resulting scatter in 
recovered STq values was consistent with the statistical 
error inferred from the SNR. Thus we infer that, in the 
ideal case of /3-model cluster profiles, uncertainty in scale 
adds 20% scatter to our uncertainty in STq. For more 
realistic cluster profiles (e.g., asymmetric), this scatter 
should be larger. 

Because the cluster counts are a falling function of 
mass (and there for of decrement amplitude), clusters de- 
tected at low significance are more likely to be lower mass 
systems with a noise-boosted decrement than a higher 
mass system with a noise-suppressed decrement. In par- 
ticular, members of the lower significance (SNR < 6) 
fraction of the ACT detections are likely to have their 
STq overestimated due to noise-boosting. In this work, 
no correction has been made for this bias. An additional 
bias arises from optimizing the detecti on significance over 
the tw o celestial coordinates and C . IVanderlinde et all 
(j201Q[ ) show that this latter bias, when referred to SNR, 
can be shown to be VSNR 2 — 3. In the case of SNR = 4, 
therefore, the optimization over location and core radius 
increases the raw decrement estimate by 10% on average. 
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Decrement measurements at 150 GHz have been pub- 
lished for a nu mber of the known clus ters in the current 
ACT sample. lHalverson etaLl (j2008l ) using the APEX 
telescope reported a decrement for th e Bullet Cluster 
-771 ± 71|aK while iPlagge et all (|2Q10f ) usi ng the South 
Pole Telescope reported —932 ± 43|liK. iPlagge et alJ 
(|2QlQh also found decrements for AS 0520 (-217^K), A 
3404 (-472(xK), AS 0592 (-529jxK), and RXC J0217.2- 
5244 (— 217|llK) all of which are in reasonable agreement 
with the 5Tq values in Table [U The central values of the 
SPT estimates are, however, systematically more nega- 
tive than the central values reported in Table [TJ 

Table [1] also contains the central Compton y parameter 
y f which is related to the central decrement STq by : 

XT' f p x -4- 1 \ 

Y^t~ B = y>0 " 4 J (! + 5 sze(x, T e )). (3) 

In this equation x — hv/ksTcMB and T e is the tem- 
perature of the electron gas. In this equation, we have 
omitted the ~ 10% contribution of the kinetic Sunyaev- 
Zel'dovich (kSZ) effect which arises due to the pe culiar 
velocity of the cluster ([Sunvaev fc Zerdovich| [l970). Be- 
cause we do not know the bulk velocity of the clusters, we 
do not account for the contribution of the kinetic effect 
which adds to the uncertainty in the derived y parame- 
ter. The y parameter itself is the product of the optical 
depth through the cluster and kBT e /m e c 2 where m e is 
the electron mass: 

Vo = -^-f / n e a T d£ (4) 
m e c 2 J 

where <jt is the Thomson cross section and isothermal- 
ity of the cluster is assumed. The relativistic correction 
Ssze(x,T c ) at 148 GHz varies from -0.03 to -0.10 for the 
ACT sample, corres ponding to cluster t emperatures in 
the range 5-15 keV (|Nozawa et al.ll2006f). We have used 
X-ray temperatures tabulated in Menanteau et~aLl ([201Qh 
as input to the relativistic correction for a subset of the 
clusters. For clusters without an X-ray temperature, we 
use a correction of -0.05 which roughly corresponds to a 
temperature of 7.5 keV, characteristic of our sample. 

For all but one case, there is no obvious correla- 
tion of the cluster decrements with strong sources. 
The one exception is Abell S0295. Approximately 
4' north-east of the cluster is ACT-S J024539-525756 
(jMarriage et al.l I2010D with a 148 GHz flux density of 
18.4 mJy (SNR=7.6). The source is also identified at 
0.84 GHz in the Sydne y University Molon glo Sky Sur- 
vey (SUMSS) catalog (jMauch et al.ll2003l ). Therefore 
the source flux density is most likely dominated by syn- 
chrotron emission. From simple geometric arguments, 
there is approximately a 10% chance that one of the the 
ACT sources would randomly fall within 4' of one of 
the twenty-three clusters, so the source-cluster proxim- 
ity may be a projection effect. Furthermore, the source's 
large angular displacement from the cluster implies that 
it is unlikely a lensed background galaxy. Finally, note 
that this source does not appear near Abell S0295 in Fig- 
ure [3] because the source was removed from the data prior 
to filtering. 

There are slight differences between the cluster prop- 
erties as presented in Table [T] and in M enanteau et al.l 
(|2010l ). The SNR and, to a lesser degree, position differ 



due primarily to t wo facts. First, the data reduction (i.e., 
the map) used for lMenanteau et al.l (|2010f ) was a prelimi- 
nary version of that used for this wor k. The map used for 
this work correspon d s to th at used in lFowler et al.l (j2Q10h 
and iMarriage et all ([2010). Second, the signal and noise 
terms used in the matched filte r for the candidate selec- 
tion in iMenanteau et all (|2010f ) are different from those 
used in this study. These differences in analysis arose 
from the need to execute the optical follow-up before the 
ACT pipeline and cluster extraction algorithm were fully 
mature. The differences in SNR generally correspond to 
approximately la fluctuations between the reductions. 

4. SIMULATED PURITY AND COMPLETENESS 

In order to assess the ratio of true detections to total 
detections (purity) and the fraction of clusters recovered 
(completeness), we run our filtering and extraction al- 
gorithm on simulations of the ACT data. The simula- 
tions were constr ucted using the millim eter signal (SZ, 
CMB, etc.) from iSehgal et all (|2010bD and noise from 
ACT difference maps. In order to avoid confusion due t o 
source-SZ correlations included in Se hgal et al.l (|2010bf ). 
the source component of those simulations was not in- 
cluded. The undetected/unmasked sources in the data 
contribute a small fraction of the rms of the filtered map 
which is dominated by wh ite detector and photon shot 
noise ([Marriage et al.ll2010D . 




SNR 



Fig. 5. — Simulated purity. Purity is defined as fraction of 
true detections in the sample. The four-pixel can didate selec- 
tion, w hich was used to produce the sample in Menanteau et al.l 
(2010), requires cluster candidates to have at least four pixels with 
SNR > 4.0. The dashed line gives the simulated purity associated 
with the four-pixel selection function. The solid line is the purity 
function corresponding to a selection function based on the single 
highest SNR cluster pixel. As expected the four-pixel purity of the 
simulated sample is in reasonable agreement with that derived from 
optica l follow-up of candidates from the data Menantea u et aTl 
(|2010D . The sample is essentially pure at SNR > 6. 

Although it is expected that, from a noise power per- 
spective, omitting the sources in this study will not sig- 
nificantly impact the analysis, there remains concern that 
millimeter sources may be correlated with clusters, either 
as cluster members or as lensed backgrounds. Studying 
573 nearby clusters with 1.4 GHz data from the NRAO 
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Fig. 6. — Cumulative completeness of the simulated sample. 
At M500C > 6 x 10 14 Mq the approximation to the total sample 
(SNR > 4.5) is ~80% complete. The completeness de grades at the 
edges of the map: taking the cleanest 230 sq-deg from lFowler et all 
( 2010), the completeness rises to 90%. 

VLA Sky Survey. iLin fc Mohrl (|2007[ ) found that the frac- 
tion of BCGs being radio-loud is higher compared to 
that of cluster galaxies of similar luminosity. Extrapolat- 
ing in frequency and redshift, this study concluded that 
roughly 10% of clusters may have 150 GHz decrements 
"filled in" by synchro tron emission from resident AGN. 
M arriage et all (j2010f ) included a cross-comparison of the 
150 GHz bright source populatio n detected by ACT and 
the 23 REFLEX X-ray clusters (jBohringer et al.l 120041 ) 
which fall within the survey area. The study cross- 
identified three clusters with 150 GHz sources. These 
three clusters were all at redshift s below z = 0.06 in con- 
trast to the higher redshift sample of this work. Lensed 
background dusty star- fo rming galaxie s may also con- 
taminate the SZ signal. 1 Wilson et al.l (|2008) found a 
bright, l ensed, dusty s ource at 1.1 mm behind the Bullet 
Cluster. IVieira et al.l ([20101 ) found that the lensed dusty 
population were characterized my maximum 150 GHz 
flux densities below 5 mJy r ange. This finding is consis- 
tent with the IWilson et all (|2008D bright Bullet Cluster 
lensed source flux of 15.9 mJy at 1.1 mm (~ 4 mJy at 150 
GHz assuming a spectral index of 2). Filtered at char- 
acteristic scale (.75') of this work, such sources would 
produce a peak brightness of 10 /iK. As the decrements 
in this work measure of order 100 /iK, contributions from 
dusty sources may affect the purity and completeness of 
faint end of the sample, but the bright end (SNR > 6) 
should be relatively immune to detection problems due 
to dusty sources. 

In order to build simulated catalogs to compare to 
that from the data, we run the filtering and detection 
pipeline described in Sec tion | 2. 21 on e a ch of s ix simula- 
tions extracted from the S ehgal et al.l (|2010bl ) sky. We 
identify decrements with clusters from the simulated cat- 
alog with M 50 oc > 2 x 10 14 M© and at redshifts z > 0.1. 
The identification radius is 3' . This radius was chosen 
by studying significant decrements detected in the sim- 
ulations and their offset from the position reported in 
the simulated catalog. By coincidence, this is also ap- 



proxi mately the size of the o ptical fields of view used 
in the iMenanteau et al.l ((2010) follow-up work. In nearly 
all cases, the arcminute scale displacements resulted from 
extended non-spherical, often double peaked, clusters. In 
these cases it was easier for noise fluctuations to move the 
detection position away from the minimum of the cluster 
potenti al as reported in the catalogs from iSehgal et al.l 

Figure [5] shows the purity as derived from the sim- 
ulations. In order to reproduc e the original selection 
function from IMenanteau et al.l (|201Q[ ) we require that 
a cluster candidate have four conjoined map pixels with 
SNR > 4.0. This requirement is intended to boost the 
purity of the sample with respect to a sample drawn with 
the requirement that only a single pixel have a high SNR. 
Figure [5] shows the purity versus SNR for both of these 
selection effects. As expected, the four-pixel requirement 
increases the purity as determined by the simulations. 
Furthermore the four-pixel pur ity agrees reasona b ly we ll 
with the purity measured in IMenanteau et all (|2010f ). 
The simulations imply that the purity of the sample with 
SNR > 6 is 95%, consistent with the rinding that all these 
clusters have optical counterparts. Furthermore, both 
the four-pixel and single-pixel selection process yield the 
same purity for the high SNR sample. Below SNR = 6, 
simulations show that the four-pixel purity falls off in a 
manner consistent with the purity determined from the 
optical follow-up. The single-pixel purity falls signifi- 
cantly faster with decreasing SNR compared to that of 
the four-pixel or optically-measured purity. We there- 
fore choose the four-pixel selection function to evaluate 
the completeness of the sample below. 

It is straightforward to simulate the cumulative com- 
pleteness as a function of cluster mass for the pure popu- 
lation at SNR > 6. As shown in Figure [6] with the dashed 
curve, we count the simulated clusters above a given mass 
and detected at SNR > 6. We then divide the number of 
detected clusters by the total number of clusters in the 
same mass range. The resulting cumulative completeness 
shows that the high SNR sample is approximately ~45% 
complete for M 5 oo c > 6 x 1O 14 M and rises to 80% com- 
plete for M500C > 1.0 x 10 15 M . One reason for incom- 
pleteness among the high mass population is that clusters 
at z < 0.3 are characterized by larger angular scales than 
their high-redshift counterparts and start blending with 
the fluctuations in the cosmic microwave background. 
Dividing the simulated sample in redshift bins of 5z = 0.2 
bins between z = 0.1 and z = 1.1, the completeness of this 
sample, in the mass range M500C > 6 x 1O 14 M , increased 
smoothly from 30% in the lowest redshift (z = 0.1-0.3) 
bin to 70% complete in the highest bin (z = 0.9-1.1). 

It is not so straightforward to reproduce the cumu- 
lative co mpleteness functio n for t he entire sample pre- 
sented in IMenanteau et al.1 ([2010) which is based on a 
pr evious reduction o f the d ata. As shown in Figure 13 
of IMenanteau et al.l (|2010f ). the follow-up observations 
were completed for a preliminary candidate list down to 
SNR = 4.2. Below SNR = 4.8, the follow-up was per- 
formed for a decreasing fraction of candidates such that, 
by SNR = 4.2, the fraction of the total candidates fol- 
lowed up was 50%. In order to approximately simulate 
the smoothly declining selection function of the follow- 
up for candidates in the range SNR = 4.2-4.8, we simply 
introduce a hard cut at SNR = 4.5 where the sampling 
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goes from 100% to 0%. 

With this prescription for modeling the entire sample, 
the simulations suggest that ~80% of all clusters with 
M^ooc > 6x 1O 14 M are detected. Figure [6] shows the full 
cumulative completeness function with a solid curve. Un- 
like in the case of the high SNR sample discussed above, 
the simulations suggest that completeness of the full sam- 
ple in the range M500C > 6 x 1O 14 M is not a strong 
function of redshift (z = 0.1-1.1). Another important 
reason for the incompleteness is that a fraction of mas- 
sive clusters are disturbed, showing elongated structure 
in the plane of the sky, and thus are not well matched to 
our f3- model filter. Conversely, the simulations suggest 
our sample likely includes high-redshift relaxed, spheri- 
cal clusters of lower mass (e.g., 4 x 1O 14 M ). Another 
feature of the data which drives the overall completeness 
down is the high level of noise towards the edge of the 
map. If the region under consideration is restricted to the 
central 2 30 square degrees u sed for the power spectrum 
study in lFowler et al.l ([2010f h then the completeness es- 
timate for clusters with M500C > 6 x 1O 14 M increases 
to -90%. 

5. SZ-MASS CORRELATION 

In this section, we perform a largely illustrative, pre- 
liminary study of the correlation of the SZ signal with 
mass for the more significant (SNR > 5) half of the 
SZ detections. For the SZ signal we use the central y' 
as described in Section [3] and listed in Table [TJ For 
the total mass, we derive M^ c using X-ray luminosi- 
ties from M enanteau et al.l (1201 Oh and the L x — M rela- 
tion from IVikhlinin et al.l (|2009f ). For all clusters but 
Abell 3128(NE), the X-ray luminosities were derived 
from the ROSAT All Sky Survey. For Abell 3128(NE), 
the available higher re solution Chandra results from 
iMenanteau et all([2010[ ) were used to help separate the 
flux from the more massive background cluster from 
the l ess massive cluster at low redshift (W erner et al.l 
120071 ). An additional scale factor of ~ 1.6 was needed 
to convert the luminositie s from the 0.1-2.4 keV band 
in IMenanteau et al.l (|2010D to the 0.5-2.0 keV ba nd re- 
quired by the correlation in IVikhlinin et al.l ([200 9) . The 
resulting M500C estimates are shown in Table [2j The 
error given for M 5 oo c only includes statistical error prop- 
agated from error on Lx- Additi onal scatter about thi s 
relation was measured to be 48% ([Vikhlinin et al.ll2009f ). 
due to evolution and non-trivial processes in the ICM as 
well as Ongoing optical and X-ray follow-up of the ACT 
cluster sample will improve these mass estimates in the 
near future. 

For self-similar evolution, the central y f is expected to 
scale as ME 2 (z) where E 2 (z) = + z) 3 + ^ A 

for a flat ACDM universe. The conventional argument 
behind this scaling is given in Appendix [A] Anticipat- 
ing such a scaling, in Figure [71 we plot y' Q values scaled 
by E~ 2 (z) against their corresponding M500C estimates. 
Values used for VLm and are 0.26 and 0.74, respec- 
tively ([Komatsu et al.lfeOlOf ). The errors plotted for y' 
are purely statistical. As discussed in Section [3j sys- 
tematic errors, primarily due to uncertainty in the core 
radius, would increase the errors on y' by at least 20%. 
Because of systematic uncertainties associated with the 
estimators for y' and M^ooc no attempt has been made 
to fit scaling relations to these data. Given this pre- 



4.0 




0.4 0.6 0.8 1.0 1.2 

M 500c (xlO 15 M ) 



Fig. 7. — SZ-Mass Correlation. Estimated Compton y parameter 
estimates (Table ^) are plotted, scaled for self-similar evolution, 
against X-ray luminosity-derived mass estimates (Table [2j) for the 
more significant half (SNR > 5) of the cluster sample. The line 
plotted is the self-similar model y' Q oc ME 2 (z) (See Appendix [All. 
The model, which is not fit to the data, agrees reasonably well with 
the ACT data. For comparison, the dashed lin e shows a simulat ion- 
motivated model with logarithmic slope 3/2 (|Motl et al] |2005l V 

caution, however, the data appear consistent with the 
self-similar model plotted as a solid line in the figure. 
An additional dashed line with y' oc M 3 l 2 E 2 (z) approx- 
imates the yo — M relation from numerical simulations 
including star fo rmation and supernova feedback from 
iMotl et al] ((2005). The data also appear consistent with 
this model. 

The SZ-Mass relation presented here is for illustrative 
purposes using central y parameter estimates and a mass 
proxy that are available now for all clusters but which 
have a large scatter. The SZ signal estimate will be im- 
proved in future analyses using additional ACT bands as 
well as additional information on clus ter gas and tem- 
perat ure profiles from the X-ray (e.g., lAndersson et al.l 
l2010h - The precise relationship between the SZ signal 
and mass needs to be determined from more precise mass 
proxies such as X-ray temperature, galaxy velocity dis- 
persion, and weak lensing. This work is in progress and 
will be reported in future publications. 

6. CONCLUSION 

ACT has identified the SZ signatures of twenty-three 
optically-confirmed galaxy clusters at a frequency of 
148 GHz in data taken during the 2008 observing sea- 
son. Of these clusters, ten are new discoveries. To- 
gether, ACT and SPT have reported twenty- five newly 
discovered systems. The Planck Satellite will soon add 
its complementary detections to the growing sample of 
SZ cluster discoveries. Due to the unique nature of the 
SZ effect, the sample is biased towards high-mass clusters 
at high redshift s. As such the sample complements cur- 
rent X-ray and optically-selected catalogues. Among the 
clusters discovered by ACT, ACT-CL J0102-4915 stands 
out as a high redshift cluster having an SZ tempera- 
ture decrement comparable to that of the most massive 
known clusters. Purity estimates from simulations of the 
cluster extraction agree with the sample p u rity d erived 
from optical follow-up in IMenanteau et al.l (|2010D . The 
simulations also estimate completeness of the sample at 
roughly 80% for M 50 o c > 6 x 1O 14 M . Over the same 
mass range, the completeness increases to greater than 
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90% if only the lowest noise half of the survey area is con- 
sidered. A comparison of the clusters' central Compton 
y parameters with X-ray luminosity-derived masses is in 
agreement with both self-similar scaling as well as scal- 
ings based on simulations which include star formation 
and supernova feedback. 

Future work with the SZ effect from galaxy clusters will 
include ACT data from multiple seasons and multiple 
millimeter bands. In addition to the southern survey 
described in this work, ACT is also surveying the celestial 
equator to complement extant and future multifrequency 
cluster work in this highly accessible strip of sky. 
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TABLE 1 

ACT 2008 148 GHz SZ-Selected Galaxy Clusters 



ACT ID RA (J2000) Dec SNR C 5T a 5T^ h y' Q z c Alternative ID 
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a The error simply reflects the SNR with no systematic uncertainty included. 

b Temperature decrement v alue with the effect of th e ACT beam deconvolved 5To x C~ 1 (6 C )- 

c Provided for reference: see Menanteau et al. ( 2010) for the sources of redshifts. 
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TABLE 2 

Cluster Mass Estimates 



ACT ID a 




b 




(xlO 14 M ) 
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a This is a subset of clusters from 
Table □ with SNR > 5. 
b Masses derived from X-ray lumi- 
no sities using the M — Lj relation 
in IVikhlinin et al.l (|2009h . The error 
only includes statistical error propa- 
gated from error on Lj. 



APPENDIX 

SELF-SIMILAR SCALING RELATION FOR Yq — M 

Assuming an isothermal ICM and taking R as the radius of the cluster, the central Compton y parameter is 
proportional to the product of the electron temperature (T), average central electron number density (no), and R: 

yo oc n TR. (Al) 

The electron density can be related through a mean molecular weight to the gas density which, assuming a standard 
gas fraction, can in turn be related to the total density (p) of dark matter plus gas. The total density is referenced to 
the background critical density (p oc p c ) which evolves with redshift as E 2 (z) = QmO- + ^) 3 + m a flat ACDM 
universe. Furthermore, the cluster radius is a function of the mass and density such that 

M\ 1/3 / M xV3 



Rk [j) k U*wJ • (A2) 

Un der the assumptions of y irialization and hydrostatic equilibrium, the temperature can be related to M and E(z) 
(e.g., iBrvan fc Normanlll998D : 

T (xM 2 ^E 2 ^(z). (A3) 
Combining these relations, the central Compton y parameter is seen to scale with mass and redshift as 

y (xME 2 (z). (A4) 



